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A comparative analysis of experimental data on electron transport in Si (100) MOSFETs in the region of high 
mobilities and strong electron-electron interaction is carried out. It is shown that electrons can be described by 
the model of a noninteracting gas with the renormalized mass and Lande factor, which allows experimentally 
verifiable predictions. 

PACS numbers: 


Although the experimental studies of the two- dimensional 
electron gas in silicon field-effect structures already have a 
history of more than half a century^, they still go on (see, 
e.g.^^). Interest in this field of research is mainly related to 
the possible metal-insulator transition^, as well as the mani¬ 
festations of strong electron-electron interaction2rdi2rJi for the 
lowest attainable electron densities on the order of 10^^ cm~^. 
In this range of electron densities, the effect of electron- 
electron interaction cannot be described in terms of a perturba¬ 
tion theory series, because the expansion parameter Vg is con¬ 
siderably larger than unity. Therefore, experimental results 
should be interpreted very cautiously and, in the first place, 
one needs to identify the common features revealed in exper¬ 
iments on samples from different manufacturers by different 
research groups using different methods of experimental data 
treatment. In this publication, a part of this program is per¬ 
formed and some further experiments are suggested. 

Almost noninteracting quasiparticles carrying a unit charge 
exist only very close to the Fermi energy. For them, one can 
define the Fermi momentum pp, determined by the number 
of free electrons introduced into the originally neutral system; 
the velocity vp', and the Lande factor g*, which characterizes 
only the quasiparticles at the Fermi surface. Each of these 
quantities can be measured experimentally. 

The most convenient and accurate experimental method to 
find the relation between the Fermi momentum and Fermi ve¬ 
locity is to determine the electron effective mass m = pp/vp 
from the analysis of the temperature dependence of the am¬ 
plitude and shape of the Shubnikov-de Haas oscillations. It is 
appropriate to begin the analysis of experimental results with 
measurements of this kind. However, it is necessary to take 
into account a number of factors that compel one to treat with 
some caution even these most reliable data, (i) The amplitude 
of quantum oscillations needs to be sufficiently small both be¬ 
cause of the limitation on the validity of the Lifshitz-Kosevich 
formula and because of the possible strong nonparabolicity of 
the spectrumi^ in the vicinity of the Fermi level. The effec¬ 
tive mass is overestimated owing to the first factor and under¬ 
estimated owing to the second, (ii) The quantum oscillation 
numbers have to be sufficiently large to satisfy the assump¬ 
tions made in the Lifshitz-Kosevich theory, (iii) The quantum 
relaxation time has to be independent of the temperature. Oth¬ 
erwise, the effective mass will be overestimated or underesti¬ 
mated, depending on the temperature ranged. 


Raw experimental data on the measurements of the effec¬ 
tive mass can be found ini^ andii. These measurements were 
carried out by different experimental groups and, furthermore, 
were performed on silicon field-effect structures made by dif¬ 
ferent manufacturers. According tc4iiS , the experimental data 
obtained for the effective mass can be well described by the 
interpolation formula 

TO = 0.205TOe(l + 0.035rs + O.OOOlbr^) (1) 

in the range 8 > > 1.5, where = 2.63(10/ns)^/^ and 

the electron density is measured in units of IO^^cto”^. 

Figure [T] shows the dependence of the effective mass on 
the electron density as a plot of nibUs/m versus Ug, where 
nib = O.IQtoe is the band effective mass. The thick solid 
line corresponds to the effective mass given by Eq.([T]), and 
points represent the experimental data fromii. In the range of 
electron densities where measurements were carried out, both 
sets of experimental data can be very well described by the 
dependence 

n, 

TO = mb- -r. (2) 

[Ug - Uc) 

Indeed, in coordinates used in EiglT] , each of these data 
sets can be fitted by a straight line with a slope of 42°. This 
corresponds to an experimentally obtained mass mb that is 
10% smaller than that known fromi. Equation @ describes 
data obtained for a large number of samples with different 
mobilities. It looks like the critical density Uc is determined 
by the interaction between electrons and depends weakly on 
the chaotic potential. 

A behavior described by Eq.(|2]i in the vicinity of the quan¬ 
tum phase transition point in a strongly interacting two- 
dimensional electron system was predicted in a number of 
publicationsi^vH. Note that Eq.(l2]i remains valid up to elec¬ 
tron densities as high as at least 10 ^^cto“^. 

The extrapolated values of the critical electron density are 
somewhat different: Ud = 0.54 * 10^^cto“^ and nc 2 = 
0.64 + lO^^CTO”^. Most probably, the 15 % discrepancy oc¬ 
curs for purely technical reasons such as different methods of 
experimental data treatment and the overestimation of the ex¬ 
perimental accuracy. Certainly, the importance of more funda¬ 
mental factors, for example, the different widths of the deple¬ 
tion layer in samples supplied by different manufacturers and. 
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FIG. 1: (Color online) Reciprocal effective mass versus the electron 
density. The vertical axis is chosen so that the corresponding depen¬ 
dence for a system of noninteracting electrons is a straight line with 
a slope of 45° passing through the origin. The thick solid line corre¬ 
sponds to the experimental results fitted by Eq. (1). Symbols and the 
straight line drawn through them represent the data from [14]. Ar¬ 
rows show the critical densities arising upon the linear extrapolation 
of the experimental data. 


FIG. 2: (Color online) Experimentally determined change in the ca¬ 
pacitance caused by the variation of the density of states from the 
value corresponding to the spin unpolarized electron system to an 
infinitely large value (symbols) and caused by the complete spin po¬ 
larization of electrons (noisy curve, B = 9.9 T). The solid horizontal 
line shows the value expected for a gas of noninteracting electrons. 
The shaded part of the plot at low electron densities corresponds to 
the region of electron localization in a magnetic field completely spin 
polarizing the electron system. 


as a consequence, the difference in the Coulomb interaction 
energy or the degree of disorder, cannot be excluded. 

Another reliably established fact is the independence of the 
effective mass of electrons in Si (100) MOSFETs on the de¬ 
gree of spin polarization of the electron system. This result 
was first obtained from the analysis of the temperature de¬ 
pendence of the Shubnikov-de Haas oscillations in the pres¬ 
ence of a magnetic field component parallel to the interfacs^i 
and recently confirmed (at least, to the first approximation) 
by independent experiments^ . This conclusion is also sup¬ 
ported by some of the raw experimental data presented in ear¬ 
lier publicationi^ and by calculations for a multivalley elec¬ 
tron system in the limiting case of weak interaction^^. 

The further interpretation of the experimental data requires 
additional assumptions that are, generally speaking, poorly 
justified. Let us assume that a system of strongly interacting 
particles can be described in terms of a Fermi gas of nonin¬ 
teracting quasiparticles with the same parameters as those at 
the Fermi level, which depend only on the electron density. 
Indirectly, this assumption is supported by a quantum Monte 
Carlo calculation, which indicated that the difference between 
the energies of a spin-polarized and a spin-unpolarized elec¬ 
tron system is proportional to the square of the degree of po¬ 
larization, similarly to the case of free electrons^. In addition, 
let us extrapolate the linear dependence in FiglUto the inter¬ 
section with the horizontal axis and further along this axis to 
zero. 

For Us > ric, the chemical potential level with respect to 
the bottom of the band will be equal to 

= UsT^ti?/2m = (us — nc)Trh'^/2mi, (3) 


According to this expression, the thermodynamic density 
of states dns/dfj, in this strongly interacting electron system 
coincides with that in the original noninteracting electron gas; 
dns/dfi = 2mb/TTh^. The validity of the latter relation can 
be tested experimentally. For this purpose, one needs to mea¬ 
sure the difference between the MOSFET capacitances in zero 
magnetic field and in those quantizing magnetic fields for 
which the thermodynamic density of states becomes infinite 
(the method is described in detail, e.g., in Ref/^). 

The results of such measurements are shown in Fig|2]. One 
can see that, in the absence of a magnetic field, the thermody¬ 
namic density of states does coincide within the experimental 
error with that of noninteracting electrons, and the above as¬ 
sumptions lead to conclusions that agree with the experiment. 

Now, let us consider the case of a magnetic field oriented 
parallel to the interface whose strength B is less than or equal 
to the field Bp required to attain the complete spin polariza¬ 
tion. The magnetic moment M per unit area of a two dimen¬ 
sional electron system is determined by the density of states 
for a fixed number of electrons and the difference between the 
energies of electrons with the spin oriented along and opposite 
to the field: 


M = 


Tfih 

27r/l^ Us — Tic 


(4) 


where ps is the Bohr magneton and g is the interaction renor¬ 
malized Lande factor. It is assumed here, in the context of 
experimental data testifying that the effective mass is indepen¬ 
dent of the parallel component of the magnetic field, that the 
critical density ric is also independent of the magnetic field. 
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FIG. 3: (Color online) Comparison of the experimental results from 
[23] (for magnetic fields of 1.5, 2, 3, 4, 5, 6, and 7 T) with (dashed 
line) the calculation according to Eq.@.). 


Using the relation 

— = 151 

dus dB' ^ ’ 

the following expression for a quantity that can be measured 
experimentally is obtained; 



FIG. 4: (Color online) Field corresponding to the complete spin po¬ 
larization of electrons versus their density. The thin line shows the 
result of calculations using the fitting line drawn through the exper¬ 
imental points in Fig. Q for g = lAgo', the thick line shows the 
expected field of complete spin polarization calculated for the pa¬ 
rameters at the chemical potential level [4, 10]. Squares and circles 
correspond to the measurements of the magnetoresistance saturation 
field [24] and the results from [23], respectively. The dash-dotted 
line is drawn according to the calculations of [20] for an ideal elec¬ 
tron system using the experimental Lande factor. 

the interpolation of the data from RefJi and the Lande factor 
g = 1.4po according to the formula 


dfj, 

'Wb 


{gBgfB 


mb 

2'Kh? (jls 



(6) 


It is noteworthy that Eq.(|6]contains only experimentally mea¬ 
surable quantities g and ric and no fitting parameters. 

Experimental data on the derivative ^ in different mag¬ 
netic fields were obtained in Ref.s22i2^ . In Eigl3 , data from 
Ref.— are compared with the curve calculated according to 
Eq.dUl with2i ric = nc 2 and g = 1.4po = 2.8. One can 
see that (i) in full agreement with Eq.l^D. the experimental 
curves in the region where they overlap do scale with the 
magnitude of the magnetic field; (ii) the dependence given 
by Eq.(|6]l qualitatively agrees with the experimental results; 
and (iii) there is a noticeable quantitative disagreement. Es¬ 
pecially troublesome is the discrepancy in the density range 
2=i=10^^cm“^ < Us < 4=1= 10^^cm“^, where it is clearly larger 
than the possible error and cannot be attributed to the finite¬ 
ness of the temperature or the inhomogeneity of the sample. 
Most probably, this discrepancy results from the use of the pa¬ 
rameters obtained for electrons at the chemical potential level 
to describe electrons within the Eermi distribution. 

Next, let us discuss the measurements of the field of com¬ 
plete spin polarization Bp for a magnetic field oriented paral¬ 
lel to the interface. The corresponding experimental data are 
shown in Eig.®. Eor comparison, the results of independent 
measurements are shown by lines. The thin solid line pass¬ 
ing slightly above the experimental points is obtained using 


g^sBp — 


risTr/i^ 

mg 


(7) 


which corresponds to the model of free electrons with the 
renormalized mass and Lande factor. The thick solid line is 
obtained from independent experiments^^, where the prod¬ 
uct gm is described by the expression 

gm = gomtil + 0.21r, -F O.OOSrf -F 0.0000045r-®). (8) 


One can see that both the experimental points and the two 
solid lines fall on straight lines with somewhat differing slopes 
at electron densities above 10^^cm“^. The difference in the 
slopes does not exceed the systematic experimental errors, and 
the extrapolation of both fitting lines to zero value of the field 
Bp yields the same critical density Uc = 0.63 * 

This value is lower than that obtained by the extrapolation of 
the straight line drawn through the experimental points, which 
is probably caused by the abovementioned overestimation of 
the energy required for the complete spin polarization of the 
electron system. 

The quantum Monte Carlo calculation of the complete spin 
polarization field^li also yields a curve close to a straight line 
in the region of interest on the {ns,Bp) plane (see EiglDi. 
However, its slope is noticeably smaller than the one obtained 
experimentally. According to Ref— , the dependence of the 
complete spin polarization field Bp on the electron density Ug 
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in an electron system with disorder is qualitatively different: 
Bp becomes zero for a finite value of Ug- This behavior re¬ 
sembles that observed experimentally. Nevertheless, even in 
this case, the calculation noticeably underestimates the energy 
required for complete spin polarization for electron densities 
exceeding 2 * . However, if the Lande factor de¬ 

termined for electrons at the chemical potential level does not 
describe electrons within the Fermi distribution and is consid¬ 
ered as an additional fitting parameter, a calculation^ using 
these two parameters (the degree of disorder and the Lande 
factor g) can fit the experimental data in FigH] quite satisfac¬ 
torily. 

The transition of an electron system subjected to a field Bp 
to a completely spin polarized state must be accompanied by 
an abrupt change in the thermodynamic density of states and, 
consequently, in the capacitance of the MOSFET structure. 
An example of a corresponding experimental plot is shown in 

FigEl 

Using Eq.®, it is easy to calculate the chemical potential 
of a partially polarized electron system in a magnetic field 
B < Bp-. 


g.{B,ns) = /i(0,n„) -f ( 


fisgB 2 mb _ 

2 irh? (ug 


ric 

- ncY ’ 


(9) 


where p(0,n„) is given by Eq.®. Then, the inverse of the 
thermodynamic density of states equals 


9/i 7rfi^ [gsgSYmbric 

dng 2mb 2'Kh'^(ng — ric)^ 


( 10 ) 


It varies most pronouncedly in the vicinity of the transition to 
the spin polarized state: 

du, , ^ , rin 

oug 2mb {ng{Bp)-nc) 

The second term in the braces in Eq. (dUl is not small, 
especially for Bp 0. Nevertheless, no experimental mani¬ 
festations of this term have been observed. Indeed, according 
to Eq. (ttUl, one would expect the occurrence of distinctly 
nonzero values of \/S.C/C\ for densities above the observed 
jump in the capacitance. This behavior is seen neither in Eigl2] 
nor in any of the experimental curves in21. 

Such a significant discrepancy arises because the quantity g, 
used in Eas. (l5l6l9ll0ll lb is the chemical potential measured 
with respect to the free electron level, while screening (and, 
thus, the capacitance) is determined by the difference between 
the chemical potential and the energy at the bottom of the two- 
dimensional electron subband. In experiments where the shift 
of the subband bottom caused by external factors is not im¬ 
portant, all above expressions are valid. The second term in 
Eq. I® equals the shift of the subband bottom, which can be 
easily seen by considering the limiting case B ^ Bp. This 
second term has to be omitted in the subsequent formulas that 
describe screening. As a result, it should be expected that the 
absolute value of the capacitance jump in the spin polarizing 


field will coincide with the horizontal line in Eig|2] Actually, 
according to the figure, this jump proves to be 25% smaller. 

Erequently, the two-dimensional electron system in 
SiMOSFETs is described in terms of another approach 
based on the Landau theory of Eermi liquids. On one hand, 
this requires the introduction of a larger number of parameters 
(amplitudes Fg , Fg , Ff instead of m and g) to be determined 
from the experiment; on the other hand, this makes it possible 
to deal only with quasiparticles within a narrow energy range 
in the vicinity of the chemical potential level. 

Since the effective mass that determines the compressibility 
and thermodynamic density of states of an electron system 
depends on two amplitudes Fg and Ff as 


^com 


m 1 + F/ 

-= mb - 

1 + FJ "l + FJ 


( 12 ) 


the weak dependence of the capacitance jump on the electron 
density (see Eig® implies^^ that Fg ~ Ff. 

It proves that the Eermi liquid parameter most sensitive to 
the interaction is Ff. According to Eqs.® and (fT2l) . F® can 
be expressed in the investigated range of electron densities as 


Ff 


r 


2 

sc 




(13) 


and is larger than 2 for an electron density of 
which corresponds to rs ~ 8. 

It is no surprise that the functional dependences of Fg and 
Ff are the same, because both amplitudes are harmonics of 
the same function F^ , and the coincidence of their magni¬ 
tudes within the experimental error means that the dependence 
of F]^ j,, on the angle between vectors k and k' is weak. 

The parameter investigated in most detail (seel^) is Fg , be¬ 
cause it is expected that it can be determined from two differ¬ 
ent experiments, i.e., from the measured values of the Lande 
factor and from the linear temperature dependence of the elas¬ 
tic relaxation tima^^. The results obtained demonstrate within 
the actual experimental accuracy that Fq depends weakly on 
the electron density in the range from 1.9* to about 

1 * 10^^cm“^.. At lower electron densities, a considerable in¬ 
crease in |Fq I takes place, in fair agreement with the result of 
Ref,26. 

Thus, any one of the methods used for experimental data 
treatment leads to qualitatively coinciding results for Fg . One 
can hardly expect better agreement because of both limited 
experimental accuracy and limited applicability of the model 
used in Ref^ and later publications. Eor example, the theory 
does not take into account in any way that the scattering poten¬ 
tial in Si MOSEETs changes radically in the range of electron 
densities under consideration. Eor the lowest densities, it is 
the screened Coulomb potential; as the density increases, it 
changes over to interface roughness scattering. Possibly, the 
experimental results obtained in Refi^, which differ some¬ 
what from those considered above, follow from the change in 
the density corresponding to the transition between different 
scattering potentials. 

In conclusion, let us once again stress the striking fact that a 
system of strongly interacting two-dimensional electrons can 
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be described with good accuracy in terms of a gas of non¬ 
interacting electrons with the renormalized mass and Lande 
factor, and this description demonstrates its predictive power. 
Hopes for further progress in the experimental investigation of 
the two-dimensional electron system in silicon MOSFETs in 
the region of low densities and high mobilities are associated 
with precision measurements of the thermodynamic density 


of states. 
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